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ABSTRACT 

We investigate star formation along the Hubble sequence using the ISO Atlas of Spiral Galaxies. Using 
mid-infrared and far-infrared flux densities normalized by K-band flux densities as indicators of recent star 
formation, we find several trends. First, star formation activity is stronger in late-type (Sc - Scd) spirals than 
in early-type (Sa - Sab) spirals. This trend is seen both in nuclear and disk activity. These results confirm 
several previous optical studies of star formation along the Hubble sequence but conflict with the conclusions 
of most of the previous studies using IRAS data, and we discuss why this might be so. Second, star formation 
is significantly more extended in later- type spirals than in early- type spirals. We suggest that these trends in 
star formation are a result of differences in the gas content and its distribution along the Hubble sequence, 
and it is these differences that promote star formation in late-type spiral galaxies. We also search for trends 
in nuclear star formation related to the presence of a bar or nuclear activity. The nuclear star formation 
activity is not significantly different between barred and unbarred galaxies. We do find that star formation 
activity appears to be inhibited in LINERs and transition objects compared to H II galaxies. The mean star 
formation rate in the sample is 1.4 M Q yr -1 based on global far-infrared fluxes. Combining these data with 
CO data gives a mean gas consumption time of 6.4 x 10 s yr, which is ~ 5 times lower than the values found 
in other studies. Finally, we find excellent support for the Schmidt Law in the correlation between molecular 
gas masses and recent star formation in this sample of spiral galaxies. 

Subject headings: galaxies: evolution — galaxies: spiral — stars: formation 
1. Introduction 
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1.1. Scientific Background 

Spiral galaxies have been recognized as sites of on- 
going star formation since the definition of the Hubble 
sequence (Hubble 1926). However the relation between 
star formation activity and the location on the Hubble 
sequence has been a continuing controversy for the past 
twenty years. 

Early optical studies comparing star formation activ- 
ity in spiral galaxies included Stauffer (1982), Kennicutt 
& Kent (1983), Kennicutt (1983). They found a gen- 
eral trend in star formation activity with Hubble type, 
with stronger star formation in Sc spirals than in Sa spi- 
rals. However, Scoville et al. (1983), using ground-based 
10 /an observations, demonstrated that star formation 
activity was not correlated with any global galaxy prop- 
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erty, including Hubble type. This laid the groundwork 
for a debate on trends in star formation along the Hubble 
sequence that continues to the present. 

With the advent of IRAS, several new studies on this 
theme were published, including de Jong et al. (1984), 
Devereux, Bccklin, & Scovillc (1987), Pompea & Rieke 
(1989), Devereux & Young (1991), Isobe & Feigelson 
(1992), Tomita, Tomita, & Saito (1996), and Devereux 
& Hamccd (1997). Most IRAS studies found no trend in 
star formation along the Hubble sequence, with the ex- 
ceptions being Devereux, Becklin, & Scoville (1987) and 
Pompea & Rickc (1989), who found that Sa galaxies are 
generally weaker sites of star formation than Sc galaxies, 
and Devereux & Hameed (1997), who found a population 
of Sa galaxies with relatively enhanced star formation ac- 
tivity. But optical and ultraviolet studies on the problem 
also continued, with Donas et al. (1987), Deharveng et al. 
(1994), and Kennicutt, Tamblyn, & Congdon (1994) all 
finding trends in increasing star formation along the Hub- 
ble sequence from Sa to Sc galaxies. To confuse matters 
further, Young et al. (1996) found trends in Ha surface 
brightness along the Hubble sequence but they argued 
that the trend did not necessarily reflect a true trend 
in star formation activity. Infrared spectroscopy of the 
C II 158 /urn line performed with ISO by Leech et al. 
(1999) also found a trend in increasing star formation 
from early- to late-type spiral galaxies, but this survey 
had a small sample of only 19 spiral galaxies. 

Clearly, opinion is sharply divided on star formation 
activity along the Hubble sequence. Some of the reasons 
for these discrepant results may include the following. 1) 
Extinction is always a concern in optical or ultraviolet 
studies, since star formation regions are always dusty. 
2) The IRAS sweep-scan technique is not ideal for mea- 
suring the flux from faint, extended sources. For faint, 
extended galaxies the IRAS point source algorithm may 
mistake the faint outer regions of the galaxy itself for 
the background and therefore undermeasure the flux. 3) 
To normalize for galaxy mass, most IRAS studies divide 
the far-infrared fluxes by B-band fluxes. Since the B- 
band fluxes represent the fluxes from the youngest, bluest 
stars in any galaxy, this tends to cancel out any trend in 
star formation activity in the far-infrared data. 4) The 
morphological classifications used are sometimes biased 
towards star formation. In particular, the Revised Shap- 
ley Ames (RSA) Catalog (Sandage & Tammann 1987) 
uses star formation activity in its classification scheme 
(Sandage & Bedke 1994), so results on star formation 
trends found by any paper using the RSA classifications 
may merely reflect the classification criteria. 

1.2. The Design of this Study 

As a fresh approach to this question, we have com- 
pleted a new survey of spiral galaxies with the Infrared 
Space Observatory (ISO) (Kessler et al. 1996) at mid- 



infrared and far-infrared wavelengths, supplemented with 
ground-based JHK photometry. The data are presented 
in the ISO Atlas of Spiral Galaxies Bendo et al. (2002) 
(henceforth referred to as Paper 1). This study incor- 
porates several features that address the concerns men- 
tioned above. 1) We minimize uncertainties related to 
extinction by working strictly with infrared data. 2) The 
ISO data has several advantages over IRAS data. ISO 
has improved sensitivity and angular resolution that en- 
ables us to distinguish nuclear from disk emission. Fur- 
thermore, ISO operates in a point and integrate mode, 
and the background is measured separately from loca- 
tions pre-selected from IRAS cirrus maps, so the flux 
from faint extended galaxies can be measured more ac- 
curately (although the flux measurements are still limited 
by background structure) . 3) We normalize mid-infrared 
and far- infrared fluxes with K-band fluxes, which rep- 
resent more accurately the total stellar population than 
B-band fluxes. 4) We use the Third Reference Catalogue 
of Bright Galaxies (RC3; de Vaucouleurs et al. (1991)) 
for galaxy classification, since the RC3 classifications are 
strictly dependent on morphological features such as the 
bulge to disk ratio and the tightness of the spiral arms 
(de Vaucouleurs 1959). 

We first summarize the sample selection, observations, 
and data reduction, the details of which can be found 
in Paper 1. Next, we discuss our diagnostics for study- 
ing star formation activity. Then we present our results 
on trends in star formation along the Hubble sequence, 
and we compare these results to previous results and dis- 
cuss the underlying mechanisms for producing the trends. 
Next, we examine the galaxies for any possible trends in 
nuclear star formation related to the presence of a bar or 
nuclear activity. We conclude by calculating quantitative 
star formation rates and gas consunsumption times and 
by comparing these results to those in the literature. 

2. Data 

2.1. The Sample 

The galaxies in this sample are a subset of a complete, 
magnitude-limited set of galaxies selected from the RSA 
Catalog. The sample comprised galaxies with Hubble 
types between SO and Sd and with magnitudes By = 12 
or brighter; galaxies in the Virgo Cluster were excluded. 
A randomly selected, subset of these galaxies was ob- 
served by ISO based on target visibility. This produced 
a total of 77 galaxies that are representative of the range 
of Hubble types in the RSA Catalog. Additional infor- 
mation on the sample and its properties is presented in 
Paper 1. 

2.2. Observations and Data Analysis 

Detailed information on the observations, the data 
processing, and the measurements themselves are pre- 
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sented in Paper 1. In short, the 60, 100, and 180 fim. 
photometry was obtained with ISOPHOT (Lcmke et al. 
1996) and processed with PIA 8.0 (Gabriel et al. 1997), 
with additional processing to correct for PSF effects. The 
12 /im images were taken with ISOCAM (Ccsarsky et al. 
1996) using the LW10 filter, which covered the range 8- 
15 /im and were processed with CIA 3.0 (Ott et al. 1997), 
followed by full processing with CIR, an implementation 
of the algorithms of Starck et al. (1999) realized by P. 
Chanial. The K-band photometric images were mostly 
taken with QUIRC at the UH 2.2 m telescope, with ad- 
ditional data from NSFCAM and SPEX at the NASA 
IRTF. 

2.3. Star Formation Diagnostics 

We seek a diagnostic of recent star formation activity. 
Both the mid- and far-infrared fluxes represent radia- 
tion from dust grains heated by the total stellar radia- 
tion field comprised of the diffuse interstellar radiation 
field (ISRF) which is dominated by the old stellar popu- 
lation, and radiation from recent star formation activity. 
The 12 /xm data includes both PAH specral line emission 
and very small grain continuum emission, both heated by 
the radiation field. While the PAH features will remain 
relatively unchanging in intensity (or decrease slightly) 
as the radiation field increases, the continuum does in- 
crease, as demonstrated by Dale et al. (2001). The 12 fim. 
data thus includes contributions from both the old stellar 
population, through the diffuse ISRF, and dust heated by 
recent star formation. Furthermore, with the exception 
of some Sd and Sm galaxies, the mid-infrared and far- 
infrared fluxes are tightly correlated among the galaxies 
in the sample, as shown in Figure 1. In particular, the 
mean ratios of the mid-infrared flux to the far-infrared 
flux in Sa, Sb, and Sc galaxies are statistically identi- 
cal, as shown in Table 1. This correlation demonstrates 
that mid-infrared fluxes are as effective as the far-infrared 
fluxes in tracing the total stellar radiation field. 

On the other hand, the K-band flux largely traces the 
old stellar population, whose thermal emission peaks at 
~ 1.6 /xm. Normalizing the 12 /im flux densities by the 
flux densities at K therefore largely cancels out the con- 
tribution from the old stellar population to the radiation 
field, so that the 12 /zm/K-band flux density ratio is a 
good diagnostic of recent star formation activity. Simi- 
larly, the far-infrared flux densities contain contributions 
from dust heated by both the diffuse ISRF and from re- 
cent star formation activity. Dividing the far-infrared 
flux densities by the K-band flux densities tends to can- 
cel the contribution from the diffuse ISRF, thus making 
far-infrared/K flux density ratios a second diagnostic for 
recent star formation activity. Therefore, we use the mid- 
and far-infrared fluxes normalized by the K-band fluxes 
as star formation indicators. 

To verify that these indicators are useful, we compare 



these flux ratios generated from the Starburst99 model 
(Lcithercr et al. 1999). We ran two test cases with the 
Starburst99 model: an instantaneous burst of star forma- 
tion that forms a fixed mass of stars, and a continuous 
burst of star formation that forms stars at a fixed rate. 
For each case, we used a Salpeter initial mass function 
with a slope of -2.35, an upper mass cutoff of 100 M©, and 
solar metallicities. The other inputs that we used for the 
model match the defaults for Starburst99. We examine 
the diagnostic Mb i - Mjf, effectively the total bolomet- 
ric flux divided by the K-band flux, which behaves like 
our normalized mid- and far-infrared fluxes. For compar- 
ison, we also examine M-Boi - Mb, or bolometric fluxes 
divided by B-band fluxes. This is akin to the far- infrared 
fluxes normalized with B-band fluxes that are often used 
in studies that examine trends in star formation along 
the Hubble sequence. 

In Figure 2 we plot M Bo i - M x and M Bo / - M B for the 
time range between 10 6 and 10 8 yr for an instantaneous 
burst of star formation, and in Figure 3 the same data 
are plotted for the continuous star formation model. The 
colors on the left side of the graphs represent younger 
stellar populations, while the colors on the right side of 
the graphs represent older stellar populations. In both 
plots, the most notable feature is the dip in the Mb i - 
Mk indicator at 10 7 yr due to the evolution of the most 
massive stars into the red supergiants that dominate the 
total K-band emission at that time. It is evident that, 
for both scenarios, the difference between recent star for- 
mation and an evolved stellar population for Mb i - 
are 1.5 - 2.5 magnitudes, and the K-band normalization 
provides one magnitude increased discrimination com- 
pared to normalization by B-band fluxes. Since the mid- 
and far-infrared fluxes are proportional to the bolometric 
luminosity, diagnostics using these fluxes divided by K- 
band fluxes are useful for tracing recent star formation, 
and are certainly preferable to diagnostics using B-band 
fluxes. Real galaxies will of course contain star formation 
regions lying within populations of older stars, and star 
formation activity may occur in a series of instantaneous 
bursts, or be continuous but time-varying. What this 
exercise does show is that the diagnostics we have cho- 
sen are effective for discriminating recent star formation 
activity. 

3. Star Formation Along the Hubble Sequence 
3.1. Star Formation Activity 

We will examine trends in star formation along the 
Hubble sequence in three regions of these galaxies. The 
global regions were defined as the regions within the inner 
135 arcsec because it matched the aperture of ISOPHOT 
observations at 60 and 100 fim. The nuclear regions were 
defined as the inner 15 arcsec, which corresponds to a 
diameter in an integer number of pixels that is wider than 
three times the full width half-maximum of the point 
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spread function in the ISOCAM images. The disk regions 
were defined as outside the nuclear regions but within the 
global regions. To study the global star formation, we 
use the far-infrared fluxes (the 40 - 220 /im fluxes from 
Paper 1) normalized by the K-band fluxes for 63 galaxies 
where we measured both fluxes within 135 arcsec. To 
study the nuclear and disk star formation, we use the 
mid-infrared (8-15 /im) fluxes normalized by the K-band 
fluxes. Nuclear star formation was investigated in 70 
galaxies in which fluxes were detected within 15 arcsec 
at > 3 a levels, and disk star formation was explored in 55 
galaxies in which fluxes within 135 arcsec were detected 
at > 3 a. The use of fixed angular apertures means that 
the physical areas sampled in each galaxy will vary with 
the distance to the galaxy. However, the flux ratios used 
here are uncorrelated with distance, so differences in the 
physical area sampled should not signifantly affect the 
analysis. 

Figures 4-6 show histograms of the star formation 
indicators for each Hubble type. A clear trend can be 
seen in these histograms, with star formation activity 
increasing from early-type to late-type galaxies. This 
trend is seen separately in global, nuclear, and disk star 
formation, although it is not as well defined in the nuclear 
regions as it is in the disk regions. 

Table 2 presents Gaussian statistics that compare the 
star formation indicators of each Hubble type. Clearly 
these statistics show quantitatively a large difference in 
the star formation activity between Sa - Sab and Sc - Scd 
galaxies. The difference between mean values of the star 
formation diagnostics for Sa - Sab and Sc - Scd galaxies 
is greater than 3 a for all regions, and the difference is 
nearly 6 a for the disk regions. Because the data may not 
necessarily obey Gaussian statistics, we also used non- 
parametric statistical methods to compare Sa - Sab, Sb - 
Sbc, and Sc - Scd galaxies. Table 3 contains results from 
applying the Kolmogorov-Smirnov (K-S) test to the di- 
agnostics for two different morphological subsets. The 
probabilities reported indicate the likelihood that the two 
populations come from the same distribution (i.e. they 
have similar values). The K-S tests show at the 98 % con- 
fidence level that the Sa - Sab and Sc - Scd nuclear star 
formation diagnostics have distinct values and at better 
than the 99 % confidence level that the Sa - Sab and 
Sc - Scd global and disk star formation diagnostics have 
distinct values. 

It has been suggested that use of RSA Catalog clas- 
sifications for studies such as these can produce spuri- 
ous correlations, since star formation is included in the 
classification scheme. The analysis above has used RC3 
classifications, but we have repeated the analysis using 
RSA classifications. We found no changes in the trends 
or their statistical significance for nuclear or global star 
formation, but the correlation of disk activity with mor- 
phological type went from ~ 6 a to ~ 8 a. Evidently the 
star formation features do make a measurable difference 



in the RSA classifications. 

These data all clearly show that star formation activ- 
ity increases along the Hubble sequence from early- to 
late-type galaxies, with a strong distinction between the 
star formation activity of Sa - Sab and Sc - Scd galax- 
ies. The trend is clear for the global, nuclear, and disk 
regions, although it is not as well defined in the nuclear 
regions. 

3.2. Spatial Extent of Star Formation 

In Paper 1 we noted that late-type galaxies usually 
contained well defined star formation structures in their 
disks while early- type galaxies often did not. Further- 
more, some late-type galaxies hosted their strongest sites 
of star formation in their disks, whereas all early- type 
galaxies produced their strongest mid-infrared emission 
from within or near their nuclei. Our qualitative conclu- 
sion was that star formation is more centralized in early- 
type galaxies and more extended in late-type galaxies. 

We now attempt to show this trend quantitatively us- 
ing the 12 fim. mid-infrared fluxes. To indicate the central 
concentration of the star formation activity, we calculate 
the logarithm of the flux within 15 arcsec divided by the 
flux within 135 arcsec, or log ( F mir(15) y Thig 

quan- 

-Fmj/?.(135 ) 

tity will be larger for a centrally concentrated galaxy and 
smaller for a more extended galaxy. (Taking the log- 
arithm of this quantity is appropriate because we are 
looking for multiplicative differences in the ratio, which 
will appear linear in logarithm space.) We perform this 
analysis for only for those galaxies where the 12 /im flux 
within the 135 arcsec aperture was detected at the 3 a 
level. Because this quantity may also be related to the 
angular width of the galaxy examined, we only use galax- 
ies with similar angular sizes, those with log (D25) values 
(from RC3) between 1.5 and 1.9. These constraints pared 
the sample down to 44 galaxies. 

Table 4 presents Gaussian statistics on the log 
( Fmir ( 15 J ) values, and Figure 7 presents histograms 

F M ir(135 ) 

showing the distribution of values for each Hubble type. 
Since this data is surely non-Gaussian, we also present 
probabilities resulting from the K-S test in Table 5. 
These data show a continuous trend from Sa to Scd 
with a 5 a significance. Alternatively, the K-S test shows 
at the 99% level that the Sa-Sab and the Sc-Scd galaxies 
differ in the spatial extent of their mid-infrared fluxes. 

3.3. Discussion 

Perhaps ironically, the ISO Atlas results are consistent 
with nearly all previous optical and ultraviolet studies - 
including Stauffer (1982), Kennicutt & Kent (1983), 
Kennicutt (1983), Donas et al. (1987), Kennicutt, Tam- 
blyn, & Congdon (1994), and Deharveng et al. (1994) - 
but inconsistent with nearly all the earlier infrared pho- 
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tometric studies. It is intriguing that our infrared ap- 
proach has confirmed the trends with Hubble type found 
in optical studies. Does this mean that the issues of 
extinction mentioned in the Introduction are actually ir- 
relevant? We suggest that extinction is an issue, and the 
infrared approach provides a statistically stronger conclu- 
sion. Unfortunately it is not possible to make quantita- 
tive comparisons because none of the optical studies that 
reported a trend with Hubble type offered any quantita- 
tive assessment of confidence. Most simply showed the 
global trend graphically as we have done in Figures 4 - 
6. 

The far-infrared spectroscopic results from Leech et 
al. (1999) are consistent with the ISO Atlas results, but 
the only two far-infrared photometric studies with results 
consistent with these are Devereux, Becklin, & Scoville 
(1987) and Pompea & Rieke (1989). Both these studies 
used RSA morphological classifications, and as we have 
argued above, trends found in samples classified with the 
RSA Catalog will at least partly reflect the classifica- 
tion criteria rather than true physical trends. The re- 
maining photometric infrared studies, including Scoville 
et al. (1983), de Jong et al. (1984), Devereux & Young 
(1991), Isobe & Feigelson (1992), Tomita, Tomita, & 
Saito (1996), and Devereux & Hameed (1997) found 
no significant correlation between star formation activ- 
ity and Hubble type. We suggest that the advantages the 
ISO measurements provide, both the angular resolution 
and the point-and-integrate observing mode, compared 
to IRAS, and the normalization by K-band fluxes rather 
than by B-band fluxes, are responsible for the different 
conclusions reached in this study compared to most of 
the previous infrared studies. 

These trends in star formation must reflect physical 
differences between early- and late- type spiral galaxies. 
These physical differences can be interpreted in the con- 
text of the results from Kennicutt (1989), which states 
that star formation is related to gas density by a sim- 
ple power law and that the gas density must exceed a 
threshold density before star formation is triggered. This 
must mean that the gas density is more likely to exceed 
the critical threshold in late-type spiral galaxies than 
in early-type spirals. In fact, the molecular and total 
gas surface density does increase from early- to late-type 
galaxies (Young & Knezek 1989), (Roberts & Haynes 
1994), and the ratio of total gas mass to dynamical mass 
also increases from early- to late-type galaxies (Young 
& Scoville 1991). This demonstrates that a typical Sc 
galaxy is more likely than an Sa galaxy to have gas densi- 
ties high enough to trigger gravitational collapse and star 
formation. Interestingly, Young et al. (1995) also demon- 
strate that molecular gas is distributed over a larger area 
in late-type spiral galaxies than in early-type spiral galax- 
ies. This suggests that our observations of more spatially 
extended star formation in late-type galaxies compared 
to early-type galaxies may be the result of the spatial 



distribution of molecular gas in spiral galaxies. 

4. Dependence of Star Formation on Other Fac- 
tors 

4.1. Relation of Nuclear Star Formation to Bars 

4-1.1. Background 

Unlike research on variation in star formation along 
the Hubble sequence, research on how bars may enhance 
star formation is reaching a consensus, but the paradigm 
is more complicated than for star formation along the 
Hubble sequence. 

Very early studies on barred galaxies, such as Sersic & 
Pastoriza (1965) and Sersic & Pastoriza (1967), demon- 
strated that barred galaxies are more likely to host en- 
hanced nuclear star formation. Studies using IRAS data, 
including de Jong et al. (1984) and Hawarden et al. 
(1986), demonstrated that barred galaxies have enhanced 
star formation compared to unbarred galaxies. These re- 
sults have been confirmed by Devereux (1987), Huang 
et al. (1996), Ho, Filippenko, & Sargent (1997a), and 
Rouusel et al. (2001), who find the enhancement only 
in early- type spiral galaxies. These results are not con- 
firmed by Isobe & Feigelson (1992), but their study did 
not subdivide the data into early- and late-type spiral 
galaxies, so they may not have been able to find the ef- 
fect seen by others. 

Although these studies tend to agree, most of them 
rely on IRAS data to detect nuclear enhancements in 
star formation, even though the IRAS detectors may not 
adequately isolate nuclear infrared emission from disk 
emission. Because the data in the ISO Atlas has several 
advantages in studying star formation, described earlier, 
we attempt to examine the issue of whether bars enhance 
star formation. Unfortunately, this sample contains very 
few early-type galaxies, where the effect is expected to 
be seen, so we cannot expect to produce a definitive con- 
clusion on this subject. 

4-1.2. Analysis 

In this analysis, we compare SA, SAB, and SB galax- 
ies for three subsets of galaxies in the sample: Sa - Sab 
galaxies, Sb - Sbc galaxies, and Sc - Scd galaxies. Since 
bars specifically affect nuclear star formation, we use 
the nuclear star formation indicator log( F ^ H ). Table 6 

presents Gaussian statistics on the log( F p^ H ) values. We 
also present probabilities from applying the K-S test in 
Table 7. 

The Gaussian statistics show very little statistical dif- 
ference between the SA, SAB, and SB galaxies for any 
subset. For Sa - Sab galaxies, we have very few galaxies 
for making a comparison, so more data would be needed 
to properly determine the presence or absence of a trend. 
For Sb - Sbc and Sc - Scd galaxies, the number of galax- 
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ies in the comparison is large enough that the statistics 
can be treated as meaningful. No difference is found be- 
tween the SA and SB types for Sb - Sbc galaxies. The 
data do show an increase in star formation from SA to 
SB galaxies for Sc - Scd galaxies, but the trend is a sta- 
tistically weak 2 a trend. The application of the K-S test 
to the data for all galaxy subsets confirms that there is 
no significant statistical difference between SA and SB 
galaxies. 

4-1.3. Discussion 

The current paradigm suggests that only in some 
barred galaxies will enhanced star formation be found, 
and then only in early-type spiral galaxies. The ISO At- 
las sample contains too few SA and SAB early-type spiral 
galaxies to make any definitive statements on how bars 
may affect nuclear star formation. Some late-type spiral 
galaxies apparently host enhanced star formation, but 
the trend is statistically weak, so we cannot make any 
definitive statements. Clearly a larger sample of galax- 
ies is required to make stronger statements on the effects 
of bars. We can however express agreement with most 
other recent results that demonstrate little enhancement 
of nuclear star formation in barred late-type galaxies. 

4.2. Relation of Nuclear Star Formation Activ- 
ity to AGN Activity 

4-2.1. Background 

The association of starbursts with active galactic nu- 
clei (AGNs) has been postulated for both Seyferts and 
low ionization nuclear emission regions (LINERs). The 
issues for these two classes of galaxies are different, how- 
ever. 

For Seyferts, the issue is how enhanced star formation 
activity is connected to the AGN activity. A few stud- 
ies, such as Rodriguez Espinosa, Rudy, & Jones (1987), 
Dultzin-Hacyan, Moles, & Masegosa (1988), Heckman et 
al. (1989), Gonzalez-Delgado & Perez (1993), Colina ct 
al. (1997), Gonzalez Delgado et al. (1998), and Oliva et 
al. (1999) have found links between the Seyfert activity 
and star formation activity, although most suggest that 
the enhancement in star formation is seen only in Seyfert 
2 galaxies. Since the ISO Atlas has superior spatial res- 
olution at mid-infrared wavelengths compared to most 
other previous infrared studies of the Seyfert/starburst 
connection, we should be able to examine whether en- 
hanced star formation activity is found in Seyferts. Un- 
fortunately, using the Ho, Filippenko, & Sargent (1997b) 
classifications for AGN, the sample contains only three 
Seyfert galaxies. 

For LINERs, the issue is whether or not LINER op- 
tical line emission is correlated with starburst activity. 
Three types of excitation mechanisms have been pro- 
posed for LINER emission: photoionization from a clus- 



ter of hot young stars (Terlevich & Melnick 1985; Shields 
1992), photoionization from a power-law source such as 
an AGN (Ho, Filippenko, & Sargent 1993), and shock 
excitation (Heckman 1980). It is not clear that only 
one excitation mechanism is responsible for producing 
all observed LINER activity as recent studies have sug- 
gested that the LINER class of galaxies may consist of 
galaxies with different central engines (Larkin et al. 1998) 
(Alonso-Herrero et al. 2000). If LINER emission is pow- 
ered by photoionization by hot young stars from a star- 
burst, the nuclei should produce relatively strong mid- 
infrared emission. If LINER emission is powered by other 
mechanisms, however, the LINER nuclei may not neces- 
sarily be sources of strong mid-infrared emission. 

To address these issues, we compare the nuclear star 
formation activity in Seyferts, LINERs, H II galaxies, 
and transition objects in our sample that have also been 
surveyed by Ho, Filippenko, & Sargent (1997b), who 
classified the galaxies' nuclear activity based strictly on 
optical spectral line ratios. (NGC 6503 is classified as 
both a transition object and a Seyfert, but we treat it 
as a transition object, which is, according to Ho, Filip- 
penko, & Sargent (1997b), the more probable designa- 
tion for the nuclear activity.) By using this one source of 
AGN classification, we ensure that the classifications are 
consistent among the galaxies in the comparison. 

4.2.2. Analysis 

We compare the various classes of galaxies using the 
nuclear star formation indicator \og{ F p' K R ), where the 
fluxes are measured for the inner 15 arcsec. Table 8 
presents Gaussian statistics on the \og( F p^ R ) values, and 
Figure 8 presents histograms showing the distribution of 
values for each Hubble type. We also present probabili- 
ties from applying the K-S test in Table 9. 

It is very apparent from the data that the Seyferts 
and H II galaxies are sites of stronger mid-infrared activ- 
ity than the LINERs and transition objects. The LIN- 
ERs and transition objects also seem to be very similar 
to each other in terms of nuclear mid-infrared activity, 
which may indicate that similar physical processes are re- 
sponsible for the generation of the spectral line emission 
in both these classes of galaxies. The Gaussian statis- 
tics show that the H II galaxies and Seyferts have higher 
log( f mim ) values than the LINERs and transition ob- 
jects. The statistical difference between the H II galaxies 
and the LINERS and transition objects is ~ 6 a. The 
results from the K-S test show that the H II galaxies have 
a distinctly different distribution of values than the LIN- 
ERs or transition objects at the 99 % confidence level. 

4.2.3. Discussion 

These results demonstrate very conclusively that the 
LINERs and transition objects in this sample are, at the 
very least, not classic starburst galaxies. They could be 
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either quiescent AGN, or galaxies with shock-dominated 
nuclear regions. The contrast in mid-infrared emission 
between the more active Seyferts and the quiescent LIN- 
ERs is also very striking, although the number of Seyferts 
in this sample is small. The data may imply that the 
presence of either circumnuclear material or circumnu- 
clear star formation will boost the AGN activity from the 
more quiescent LINER level to the more active Seyfert 
level. However, more work is needed to investigate the 
connection, specifically a more careful study of the cir- 
cumnuclear regions of Seyferts and LINERs as well as a 
mid-infrared survey of all Seyferts and LINERs from the 
sample in Ho, Filippenko, & Sargent (1997b). 

It is noteworthy, but perhaps not surprising, that the 
H II galaxies and Seyferts are indistinguishable using 
these photometric diagnostics. This is, after all, the un- 
derlying cause of the "starbursts- monsters" controversy 
(Heckman et al. 1983). 

5. Star Formation Rates and Molecular Gas 
Content 

5.1. Global Star Formation Rates 

To measure the range of star formation activity in 
the ISO Atlas sample, we have converted the far-infrared 
fluxes into quantitative star formation rates by assuming 
that all the far-infrared flux comes from star formation 
regions. We use the formula 

SFRiMeyr- 1 ) = 1.76 x IQ- W L FIR (L & ) (1) 

where Lfir is the 8 - 1000 ^m flux (Kennicutt 1998b). 
Both Sauvage & Thuan (1992) and Kennicutt (1998a) 
have emphasized that direct conversion of far-infrared 
fluxes to star formation rates relies on the assumptions 
that far-infrared fluxes are emitted from dust heated 
by star formation regions and that the dust is optically 
thick. Dust heated by an interstellar radiation field gen- 
erated by old stars will also produce far-infrared emis- 
sion, leading to star formation rates that are calculated 
to be higher than they really are. Nonetheless, we will 
use this conversion factor as an estimate of the global 
star formation rates in these galaxies. 

First we converted the global 40 - 220 /im fluxes 
(i.e. the fluxes from within the inner 135 arcsec) to 
8 - 1000 yum fluxes using the conversion factor 1.192 ± 
0.235. Following a method similar to the one described 
in Stickel et al. (2000), this conversion factor and its un- 
certanties were calculated by assuming that the emission 
in the 8 - 1000 /im range is dominated by a single black- 
body component with temperatures ranging between 20 - 
40 K and emissivities with wavelength dependence rang- 
ing from A to A~ 2 . Then luminosities were calculated, 
and the above formula was applied to convert the lu- 
minosities into star formation rates. Figure 9 shows the 
range of star formation rates. The average star formation 



rate is 1.4 M Q yr -1 , with a maximum of 8.5 M yr -1 and 
a minimum of 0.0035 M Q yr -1 . These results are affected 
by Malmquist bias because the sample was magnitude- 
limited. Also, the validity of the conversion factor has 
not been tested through observations. Nonetheless, these 
star formation rates are consistent with the spiral disk 
star formation rates listed in Kennicutt (1998a). 

5.2. Gas Consumption Times 

Since we have estimated global star formation rates 
for these galaxies we should also examine whether the 
galaxies contain sufficient gas to maintain this level of 
star formation activity over a Hubble time. This issue 
has been investigated previously by Larson, Tinsley, & 
Caldwell (1980), Kennicutt (1983), Sandage (1986), 
Donas et al. (1987), Kennicutt, Tamblyn, & Congdon 
(1994), and Devcreux & Hameed (1997). All of these 
studies have found gas consumption times of a few times 
10 9 yr, which is shorter than a Hubble time. They usu- 
ally suggest that gas recycling is responsible for increas- 
ing the gas consumption time, although only Kennicutt, 
Tamblyn, & Congdon (1994) have modeled this effect. 
However, most of the above studies have relied on optical 
or ultraviolet data for determining star formation rates. 
Dust extinction in these wavebands may lead to calcu- 
lated star formation rates that are artificially low, which 
then may lead to calculated gas consumption times that 
are artificially high. Wc attempt to address this issue 
with our far-infrared data, despite the concerns about us- 
ing the data to calculate star formation rates mentioned 
above. 

To calculate gas consumption times, we will use the 
"Roberts time" 



suggested by Kennicutt, Tamblyn, & Congdon (1994). 
For gas masses, we used the CO data from the FCRAO 
Extragalactic CO Survey (Young et al. 1995) to estimate 
the gas masses within 135 arcsec apertures. We calcu- 
late gas consumption times for the 20 galaxies in the ISO 
sample for which the CO intensities were detected at the 
3 a level in the central beam. The average gas consump- 
tion time we calculate is 6 x 10 8 yr, with a maximum of 
1.3 x 10 9 yr and a minimum of 3 x 10 8 yr. These gas 
consumption times seem to be independent of Hubble 
type or the star formation metrics we used above. 

The gas consumption times we calculate are short 
compared to the few times 10 9 yr values found in most 
other studies. It is possible that the high star formation 
rates we calculate are more accurate than those calcu- 
lated using data from optical and ultraviolet wavebands 
because of dust extinction issues. Alternatively, our star 
formation rates could plausibly be high for several rea- 
sons. There may be contributions to the far-infrared flux 
from dust heated by the old stellar population that could 
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increase our calculated star formation rates. Also, if we 
had used a different conversion factor to convert infrared 
flux to star formation, for example that given by Scoville 
& Young (1983), we would have found gas consumption 
times on the order of 10 9 yr, which would be more in line 
with most previous results (most closely with Devereux 
& Hameed (1997), who also calculated gas consumption 
times from far- infrared data). However, our global star 
formation rates are comparable to those listed in Kenni- 
cutt's review (Kennicutt 1998a), so our star formation 
rates may be correct. In summary, multiple factors may 
be responsible for the slightly low gas consumption times 
we find compared to most previous results, but because 
of uncertainties in the assumptions, we cannot argue that 
our results are more accurate than most previous calcu- 
lations. 

5.3. Testing the Schmidt law 

In discussing the trend in star formation with Hubble 
type, we suggested that it resulted from the difference 
in molecular gas content between early- and late-type 
galaxies. If so, for galaxies where we have both infrared 
and CO data, we should be able to find the power law re- 
lation between star formation and gas density discussed 
by Schmidt (1959) and Kennicutt (1989). We investi- 
gate this with the molecular gas mass data from Young 
et al. (1995) and do the comparison with star formation 
activity for global, nuclear, and disk regions separately. 
For the global analysis we use the star formation rates 
found above, and for the nuclear and disk analysis we 
use the 12 and 60 /im flux densities, which have spatial 
resolutions that are equivalent to or better than the CO 
data. Because we were limited by the resoulution of the 
CO data, we used 45 arcsec apertures to represent the 
nuclear regions (instead of the 15 arcsec aperutres used 
in the analysis of star formation trends). 

Figure 10 shows the gas masses versus star formation 
rates normalized by the square of the distance to the 
galaxies. We use this normalization because both quan- 
tities are effectively derived from luminosities; otherwise 
there may be a spurious correlation since one is effec- 
tively plotting D 2 vs. D 2 . The figure shows that the two 
quantities are very well correlated and can be represented 
by a simple power law. 

Figures 11 and 12 show the normalized gas masses 
versus 12 and 60 /im flux densities for the nuclear regions 
of the galaxies. Most of the galaxies lie very nicely along 
this trend. The exceptions are NGC 1569 (the outlier 
near the upper left of the figures) and NGC 5236 (the 
outlier near the upper right of the figures). NGC 1569 is a 
nearby dwarf galaxy with strong star formation activity, 
so it may not necessarily follow the same trend as most 
spiral galaxies. 

Figures 13 and 14 show the normalized gas masses 
versus 12 and 60 (im flux densities for the disk regions 



between the 45 and 135 arcsec diameters of galaxies. The 
trends are not as well defined as for the nuclear regions, 
particularly at 12 /im. The scatter, particularly at low 
gas masses, may represent the failure of gas densities to 
exceed the threshold necessary for star formation as dis- 
cussed by Kennicutt (1989). 

In summary, these comparisons demonstrate an excel- 
lent Schmidt law relating molecular gas and cither star 
formation rates or far-infrared flux densities. The trend 
is seen globally as well as in separate nuclear and disk 
regions, with occasional deviations from the trend found 
in unusual objects or in objects where gas masses fail to 
exceed the threshold for star formation in their disks. 

6. Summary 

Using the ISO Atlas data for an optically-selected, 
magnitude-limited sample of spiral galaxies, combined 
with normalization by K-band flux densities, we have ex- 
amined the evidence for star formation in three regions 
of these galaxies: integrated over the entire 135 arcsec 
of the galaxies (global), the central 15 arcsec (nuclear), 
and in the disk itself excluding the nuclear region. We 
have found strong trends in star formation activity along 
the Hubble sequence in all three regions with statistical 
significance of ~ 3-6 a. We also show that the distribu- 
tion of star formation regions is more spatially-extended 
in late-type galaxies than early-type galaxies. 

These results confirm most previous optical and ultra- 
violet investigations of star formation along the Hubble 
sequence, but do not confirm the findings of an equally 
large number of such studies using infrared (mostly 
IRAS) data, which have found no correlations with Hub- 
ble type. We suggest that the superior angular resolu- 
tion in the ISO Atlas data, combined with use of K-band 
flux densities for normalization, has resolved this long- 
standing discrepancy between optical and infrared stud- 
ies. 

We argue that the differences in star formation activ- 
ity along the Hubble sequence may be understood pri- 
marily in terms of the differences in the molecular gas 
density and distribution. We confirm that star forma- 
tion is well-correlated with molecular gas mass. We sug- 
gest that the strength of the star formation activity in- 
creases as the gas density increases along the Hubble se- 
quence. Similarly, we suggest that the spatial extent of 
star formation activity increases as the spatial extent of 
the molecular gas increases. 

We found no enhancement of nuclear star formation 
in barred galaxies compared to unbarred galaxies, in con- 
trast with many previous results. This may be because 
the effect is seen chiefly in early-type galaxies and we 
have few early-type SA and SAB galaxies in the sample. 
We do find that star formation seems to be inhibited in 
LINERs and transition objects compared to H II galax- 
ies, suggesting that LINER emission is not powered by 



8 



recent starburst activity (at least in nearby galaxies). 

From the far-infrared fluxes we estimate star forma- 
tion rates of 1.4 M Q yr -1 , which is in agreement with pre- 
vious studies of similar spiral galaxy samples. However, 
we find gas consumption times of 6 x 10~ 8 yr, which is 
about a factor of 5 lower than most previous studies. We 
do find excellent correlation between molecular gas den- 
sities and recent star formation activity in accord with 
the Schmidt Law. 

The ISOPHOT data presented in this paper were re- 
duced using PIA, which is a joint development by the 
ESA Astrophysics Division and the ISOPHOT Consor- 
tium led by MPIA, with the collaboration of the Infrared 
Processing and Analysis Center (IPAC). Contributing 
ISOPHOT Consortium institutes are DIAS, RAL, AIP, 
MPIK, and MPIA. 

The ISOCAM data presented in this paper were anal- 
ysed using 'CIA', a joint development by the ESA As- 
trophysics Division and the ISOCAM Consortium. The 
ISOCAM Consortium is led by the ISOCAM PI, C. Ce- 
sarsky. 

This research has been supported by NASA grants 
NAG 5-3370 and JPL 961566. 
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Fig. 1. — Plot of the mid-infrared fluxes within 135" 
versus the far-infrared fluxes within 135" for all galaxies 
detected in both wavebands. The line shows the best 
fitting power law. 
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Fig. 3. — Plots of two star formation indicators following 
the onset of continuous star formation, as modeled by 
Starburst99. The solid line represents Mb i - M^- and 
the dashed line represents Mb i - Mg. 



Fig. 2. — Plots of two star formation indicators following 
an instantaneous burst of star formation, as modeled by 
Starburst99. The solid line represents Mb i - Mk and 
the dashed line represents Mb i - Mg. 
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Fig. 6. — Histograms of the disk star formation diagnos- 
tic log(-^Jrp) for subsets of the galaxies in the sample 
binned according to Hubble type. 
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Fig. 7. — Histograms of the diagnostic of the central con- 
centration of star formation, log ( Fmir ( 15 J \ f or subsets 

Fmir(135 ) 

of the galaxies in the sample binned according to Hubble 
type. Only galaxies where 1.5 < log(D 2 5) < 1.9 were used 
in these histograms. The diagnostic is high for galaxies 
that are centrally concentrated and low for galaxies that 
are extended. 
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Fig. 8. — Histograms of the nuclear star formation diag- 
nostic \og( F p^ R ) for subsets of the galaxies in the sample 
based on the nuclear activity classifications in Ho, Filip- 
penko, & Sargent (1997b). Note the relatively lower 
star formation activity in LINERs and transition objects 
compared to H II galaxies and Seyferts. 
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Fig. 9. — Histogram of the distribution of star formation 
rates for all galaxies in the sample where far-infrared lu- 
minosities were calculated. 
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Fig. 10. — Plot of the global star formation rate ver- 
sus molecular gas mass for galaxies in both the ISO and 
FCRAO samples. Both quantities have been normalized 
by the distance squared for reasons explained in the text. 
The line shows the best fitting power law. 



15 



10000 



c 
o 

o 
E 

CN 



V) 

c 

01 
Q 

x 

13 



1000 




100 



M gos / D 2 (Solar masses / Mpc 2 ) 



Fig. 11. — Plot of the 12 /an flux densities versus molec- 
ular gas mass for the nuclei of galaxies in both the ISO 
and FCRAO samples. The gas masses in this figure and 
Figures 12-14 have been normalized by the distance 
squared for reasons explained in the text. The lines in 
this figure and Figures 12 - 14 show the best fitting power 
laws. 
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Fig. 13. — Plot of the 12 /xm flux densities versus molec- 
ular gas mass for the disks of galaxies in both the ISO 
and FCRAO samples. 
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Fig. 12. — Plot of the 60 /<m flux densities versus molec- 
ular gas mass for the nuclei of galaxies in both the ISO 
and FCRAO samples. 
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Fig. 14. - Plot of the 60 /im flux densities versus molec- 
ular gas mass for the disks of galaxies in both the ISO 
and FCRAO samples. 
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Table 1 

Ratio of Mid- to Far-Infrared Fluxes within 135" Apertures for Different Hubble Types 



Type 


No. Galaxies 


Mean log(-S^) 

ov J^FIR ' 


E - SO/a 


6 


-1.05 ± 0.05 


Sa - Sab 


7 


-0.94 ± 0.02 


Sb - Sbc 


17 


-0.94 ± 0.03 


Sc - Scd 


21 


-0.94 ± 0.03 


Sd - Im 


6 


-1.24 ± 0.14 


All 


57 


-0.98 ± 0.02 



Table 2 

Star Formation Indicators for Different Hubble Types 



Type Global Star Formation Nuclear Star Formation Disk Star Formation 

No. Galaxies Mean logt^ 3 -) No. Galaxies Mean log(- E f ia ) No. Galaxies Mean log(- E f ia ) 



E - SO/a 


11 


-0.03 ± 0.10 


12 


-0.92 ± 0.11 


6 


-0.47 ± 0.16 


Sa - Sab 


7 


0.21 ± 0.09 


8 


-0.69 ± 0.10 


5 


-0.06 ± 0.14 


Sb - Sbc 


15 


0.69 ± 0.09 


17 


-0.33 ± 0.13 


17 


0.58 ± 0.09 


Sc - Scd 


23 


0.67 ± 0.05 


26 


-0.30 ± 0.07 


21 


0.85 ± 0.06 


Sd - Im 


7 


0.66 ± 0.11 


7 


-0.13 ± 0.14 


6 


0.87 ± 0.12 


All 


63 


0.50 ± 0.05 


70 


-0.44 ± 0.06 


55 


0.54 ± 0.07 
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Table 3 

Results from Applying the K-S Test to Star Formation Indicators 





Data Sets Used 


Global log(^a-) 


Probability a 
Nuclear logf^&f^) 


Disk log(^) 


Sa 


- Sab and Sb - Sbc 


0.027 


0.063 


0.034 


Sa 


- Sab and Sc - Scd 


0.0010 


0.021 


0.0014 


Sb 


- Sbc and Sc - Scd 


0.33 


0.52 


0.083 



a The probability of two data sets coming from the same distribution, as determined 
by applying the K-S test. 



Table 4 

Mid-Infrared Spatial Extent Indicators for Different Hubble Types 



Type 


No. Galaxies 


Mean log( * UI *^,l ) 


E - SO/a 


4 


-0.49 ± 0.08 


Sa - Sab 


5 


-0.69 ± 0.07 


Sb - Sbc 


15 


-0.99 ± 0.10 


Sc - Scd 


15 


-1.21 ± 0.06 


Sd - Im 


5 


-1.11 ± 0.14 


All 


44 


-1.00 ± 0.05 
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Table 5 

Results from Applying the K-S Test to Mid-Infrared Spatial Extent Indicators 





Data Sets Used 


Probability a 


Sa 


- Sab and Sb - Sbc 


0.23 


Sa 


- Sab and Sc - Scd 


0.011 


Sb 


- Sbc and Sc - Scd 


0.095 



a The probability of two data sets 
coming from the same distribution, as 
determined by applying the K-S test. 



Table 6 

Mid-Infrared Nuclear Star Formation Indicators for SA, SAB, and SB Galaxies 



Type No. Galaxies Mean log(^f^) 



SAa- 


SAab 


2 


-0.81 ± 0.38 


SABa- 


SABab 


1 


-0.77 


SBa- 


SBab 


5 


-0.63 ± 0.12 


SAb- 


SAbc 


6 


-0.18 ± 0.19 


SABb- 


SABbc 


7 


-0.46 ± 0.24 


SBb- 


SBbc 


4 


-0.31 ± 0.28 


SAc- 


SAcd 


11 


-0.47 ± 0.09 


SABc- 


SABcd 


9 


-0.25 ± 0.12 


SBc- 


SBcd 


6 


-0.08 ± 0.14 
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Table 7 



Results from Applying the K-S Test to Mid-Infrared Nuclear Star Formation Indicators for SA and 

SB Galaxies 





Data Sets Used 




Probability a 


SAa- 


SAab and SBa - 


SBab 


0.96 


SAb - 


SAbc and SBb - 


SBbc 


0.98 


SAc - 


SAcd and SBc - 


SBcd 


0.16 



a The probability of two data sets coming 
from the same distribution, as determined by 
applying the K-S test. 



Table 8 

Mid-Infrared Nuclear Star Formation Indicators for Galaxies with Different Nuclear Activity 

Types 



Type 


No. Galaxies 


Mean log(^^s) 


H 11 Galaxies 


30 


-0.26 ± 0.06 


Seyfert 


3 


-0.34 ± 0.25 


LINER 


8 


-0.88 ± 0.11 


Transition Objects 


10 


-0.83 ± 0.12 



Table 9 

Results from the K-S Test to Mid-Infrared Nuclear Star Formation Indicators for Galaxies with 

Different Nuclear Activity Types 



Data Sets Used 


Probability a 


H II Galaxies and Seyferts 


0.999 


H II Galaxies and LINERs 


0.0023 


H II Galaxies and Transition Objects 


0.013 


Seyferts and LINERs 


0.21 


LINERs and Transition Objects 


0.71 



a The probability of two data sets coming from the 
same distribution, as determined by applying the K-S 
test. 
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